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We compare quark stars made of color-superconducting quark matter to normal-conducting quark 
stars. We focus on the most simple color-superconducting system, a two-flavor color superconductor, 
and employ the Nambu-Jona-Lasinio (NJL) model to compute the gap parameter and the equation 
of state. By varying the strength of the four-fermion coupling of the NJL model, we study the mass 
and the radius of the quark star as a function of the value of the gap parameter. If the coupling 
constant exceeds a critical value, the gap parameter does not vanish even at zero density. For coupling 
constants below this critical value, mass and radius of a color-superconducting quark star change 
at most by ~ 20% compared to a star consisting of normal-conducting quark matter. For coupling 
constants above the critical value mass and radius may change by factors of two or more. 



I. INTRODUCTION 

At sufficiently high densities and sufficiently low temperatures quark matter is a color superconductor In nature, 
color-superconducting quark matter could exist in the interior of compact stellar objects such as neutron or quark 
stars. Among the best known properties of compact stellar objects are their masses and radii. The question then is 
whether these observable properties allow to decide if compact stellar objects contain, or are even completely made 
of, color-superconducting quark matter. To this end, one has to compute these properties for stars containing color- 
superconducting quark matter and compare them to the corresponding ones for stars containing normal-conducting 
quark matter. 

This question has recently triggered a lot of activity 0, S B IE S S Hi IS ^i Oi Qi ^| ■ Pure quark stars as well as 
hybrid stars were considered, both with two-flavor color-superconducting quark matter as well as with quark matter 
consisting of three flavors in the color-flavor-locked phase Kljj. All these investigations are based on variants of the 
NJL model where fermions interact via a four-point vertex. The coupling strength is adjusted to be in agreement 
with hadron phenomenology at zero quark-chemical potential fi. This leads to color-superconducting gap parameters 
(j) of the order of 100 MeV (l3. E3| . For reasonable values of the parameters entering the equation of state, such as the 
MIT bag constant B and the strange quark mass nis , the result of these studies is that mass and radius of a compact 
stellar object change by ~ 20%, if it contains color-superconducting instead of normal-conducting quark matter. 

In this paper, we consider a different question. We ask how large the color-superconducting gap parameter has to be 
in order to see substantial changes in mass and radius of a compact stellar object. As the transition to hadronic matter 
introduces another degree of freedom which may either mask Q or enhance the effects of color superconductivity, 
we do not consider hybrid stars, but focus exclusively on pure quark stars. We also consider the most simple color- 
superconducting state, nar aely quark matter with two flavors in the so-called 2SC phase, although this state may not 
be the most favorable one (l7l |. 

This paper is organized in the following way. In Sec. ^ we derive the gap equation and the equation of state for 
two-flavor color-superconducting quark matter using the Cornwall- Jackiw-Tomboulis (CJT) formalism While this 
formalism is equivalent to other approaches to derive the gap equation and the equation of state, it is nevertheless the 
most elegant way. Moreover, it also provides a general framework that allows one to go beyond the standard mean- 
field approximation (although this direction is not pursued in this work). In addition, it accounts for the possibility 
of non- vanishing gluon background fields generated by condensation of quark Cooper-pairs. Our derivation presented 
in Sec. ini puts special emphasis on this point, which has previously been neglected in the the derivation of the gap 
equation. In Sec. IIIII we compute the masses and radii of quark stars via the Tolman-Oppenheimer-Volkoff (TOV) 
equation. Section Hvl concludes this paper with a summary of our results. 

Our units are h — c — ks — I- The metric tensor is 5^1/ = diag (1, — 1, — 1, — 1). Four- vectors are denoted as 
K^' = (fc°,k), where k is a three- vector with modulus A; = |k| and direction k = k/fc. We work in the imaginary- 
time formalism, i.e., the space-time integration is defined as J-^ = Jp^^dr d^x, where r is Euclidean time, T 
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is the temperature, and V the three-volume of the system. Energy-momentum sums are written as T/VJ2k ~ 
T J d'^k/ (2Tr)^, where the sum runs over the Matsubara frequencies w„ = 2mrT for bosons and ujn = (2n + 1)ttT 
for fermions, respectively. 



II. EQUATION OF STATE AND GAP EQUATION 



For color-superconducting matter in the 2SC phase, the Lagrangian is given by 

1 

4 

where 



abc A b 



(1) 



(2) 



is the gluon field strength tensor; is the vector potential for the gluon field, f are the structure constants of 



SU{3)c, and g is the strong coupling constant. The fermion fields tJj are {ANcNf 
the Dirac structure we choose the following basis in color-fiavor space: 



24)-dimensional spinors. Suppressing 



(3) 



The Dirac conjugate spinor is defined as tp = 7o- The covariant derivative is given by = 9^ — igA'^^Ta, where Tq 
are the generators of SU{i)c, suitably generalized to our 6-dimensional color-fiavor basis (PJ. The quark mass matrix 
is TO = dia.g{rnu,rnd,mu,rnd,rnu,rnd)- When computing quark star properties in Sec. lIIIL we also include electrons in 
order to achieve electrical neutrality. In some cases we also add non-interacting strange quarks. 

In the treatment of superconducting systems it is advantageous to double the fermionic degrees of freedom by 
introducing Nambu-Gor'kov spinors 



(4) 



where 4'c = Cip'^ is the charge-conjugate spinor; C is the charge-conjugation matrix. In this basis, the tree- level 
action can be written as 



X J x^ 



where 











(5) 



(6) 



<c _ 



is the tree-level propagator for Nambu-Gor'kov fermions. Here we introduced the charge-conjugate covariant derivative 
d^ + igA'^Tj' . The delta function is defined as S^^^X — Y) — dijx — Ty) 5^^\-x.—y). The quark-chemical potential 
diag(^^, /xj?, Mg J Mo' Mb J Mb)- The chemical potential for quarks of color i and fiavor / can be represented 



matrix is ji 
as 



(7) 



where /ie is the electro-chemical potential, is the electric charge (in units of e) of quark fiavor /, and /is and 
/xs are the color-chemical potentials associated with the diagonal generators and of SU{^)c- While /x controls 
the quark number density, /te, M3) ^ind /xg have to be introduced to ensure electric- and color-charge neutrality. If 
the SU{?>)c color symmetry is not broken, the color-chemical potentials have to vanish, f^s = fJ,s = 0, otherwise they 
would break SU{3)c explicitly. However, when the color symmetry is broken by a color-charged quark Cooper-pair 
condensate, /is and /tg do not need to be zero. We shall come back to this issue below. 
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The eflective action in the CJT formaUsm reads 

r[^,^5j,A,S,D] = l[^,^,A]-hTh-iD-' -^Tt{D-^D-1) 

+ ^TY\nS-^ + ^TT{S^'S-l)+T2[^,'^,A,S,D] . (8) 

The quantities D and S are the full gluon and quark propagators, respectively. The inverse tree- level quark propagator 
Sq^ was introduced in Eq. ©. Correspondingly, D^^ is the inverse tree-level gluon propagator. The traces run over 
space-time, Nambu-Gor'kov, color, flavor, and Dirac indices. The factor 1/2 in front of the fermionic one-loop terms 
compensates the doubhng of the degrees of freedom in the Nambu-Gor'kov basis. The functional r2 is the sum of all 
two-particle irreducible (2PI) diagrams. It is impossible to evaluate all 2PI diagrams exactly. However, the advantage 
of the CJT effective action ijHJ is that truncating the sum r2 after a flnite number of terms still provides a well-deflned 
many-body approximation. Here we only include the sunset-type diagram shown in Fig. 



FIG. 1: The sunset-type diagram. 



r2 = - ^ / Tr [rj^ Six, Y) Tl S{Y, X)] Df,{X, Y) , (9) 

4 Jx,Y NG,cJ,s 

where the trace now runs only over Nambu-Gor'kov, color, flavor, and Dirac indices. The Nambu-Gor'kov vertices are 
defined as 

Later on, we shall approximate the gluon-exchange interaction between quarks by a point-like four-fermion coupling. 
This effectively removes dynamical gluon degrees of freedom, such that we do not need to worry about gauge fixing 
or possible ghost degrees of freedom. Therefore we already omitted the latter in Eq. (jHJ. 

The stationary points of the effective action lISl determine the expectation values of the one- and two-point functions, 



ST ST 



ST ST ST 

IA^='' SD='^ Is-'- 



(11) 



The first two equations yield the Dirac equation for the fermionic fields 5* and ^ in the presence of the gluon field 
A^. The solution is trivial, since fermionic, i.e. Grassmann- valued, fields do not have a (c-number) expectation value. 
The third equation is the Yang-Mills equation for the gluon field. 



VfF^'^iX) 



SAliX) 



\Tr{D-^D-S,'S)~T, 



(12) 



where I?"'' = duS°'^ — is the covariant derivative in the adjoint representation. The first two terms on the 

right-hand side are the contributions from gluon and fermion tadpoles j2^. The functional derivative with respect to 
acting on the trace is nontrivial because of the dependence of the inverse tree- level propagators D^^ and Sq^ on 
the gluon field, cf. Eq. ©. The last term is non-zero if T2 contains 2PI diagrams with an expHcit dependence on A^. 
It vanishes in our approximation ^ for r2. As shown in Ref. [l^l the solution of the Yang-Mills equation in the 2SC 
phase is a constant background field A° ~ g^oS""^. This background field acts like a color-chemical potential and 
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provides the color-charge neutrahty of the 2SC phase [231 • Later on, we shall remove the gluon degrees of freedom by 
approximating the non-local gluon exchange with a point-like four-fermion coupling. The constant background field 
then disappears from the treatment, and the color-chemical potential jig assumes the role of the background field 
to ensure color neutrality. 

The fourth equation Hll|l is the Dyson-Schwinger equation for the gluon propagator, 

D-'Z{X,Y)^D^ZiX,Y) + U^/,iX,Y) , (13) 



where 

n^/,iX, Y) = -2 ^j^u^^^^x^ = y Tr [r;^ S{X, Y) n S{Y, X)] (14) 



is the gluon self-energy. Since we shall approximate gluon exchange by a four-fermion coupling, we do not need to 
solve the Dyson-Schwinger equation for the gluon propagator. 

The fifth equation (jllll is the Dyson-Schwinger equation for the quark propagator, 

S-\X,Y)^Sa\X,Y) + ^{X,Y) , (15) 

where 

m, y) = 2 J^^^ = S{X, Y) Tl Df^{Y, X) (16) 



is the quark self-energy. Assuming translational invariance, in momentum space the Dyson-Schwinger equation reads 

T 
V 



S-\K) = S^\K) + nK). m)^-9'^T.^^aSiQ)nD;l{K~Q). (17) 



Let us introduce the Nambu-Gor'kov matrices 




where 



[G+] ' [K) = r{K, + gA^Ta) + Ato - m , (19a) 



[Go ] ' (K) = riK, - g^^Tj) - A7o - m (19b) 

are the inverse tree-level propagators for particles and charge-conjugate particles, respectively. The quantities in Eq. 
itTsjl are the normal self-energies for particles and charge-conjugate particles, while $^ are the so-called anomalous 
self-energies. The self-energies are related via S"(i^) = C[J:+{-K)fC-'^ and <^>-{K) = -fo[<^>+{K)]''jo. With the 
definitions (|18ll . in Nambu-Gor'kov space the Dyson-Schwinger equation Ijl7|l has the solution 



\arriT 



where 



G± = |[G±]-Vs±~<i>^([Gf]-VsT)''<i>±| \ (21a) 
5± = - ([Gjj^Vs'FV^^iG* . (21b) 



Here G± are the propagators for quasiparticles and charge-conjugate quasiparticles, respectively, while 5* are the 
so-called anomalous propagators. 

The gap equation for the color-superconducting gap parameter can be deduced from the (21)- or (12)-components 
of the Nambu-Gor'kov self-energy Ijl8ll . 

<^+[K) ^g'^J2 ^'^a S+(g) Df,{K ~ Q) , (22a) 
Q 

$-(/0 = ,92|^^A^T,S-(Q)7'^r,^Df,(i^-Q) . (22b) 
Q 
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It is sufficient to consider Eq. H22a|l . because Eq. Ij22b|l follows from the relation $ (K) = jct[^~^ (K)]^ jo . In the 
color-fiavor basis ^ the gap matrix in the 2SC phase reads 

/ 
A± 
±^ 

A^^ 

\ 

Here we have used the fact that red up-quarks form Cooper pairs with green down-quarks, with a (Dirac space) 
gap matrix Af, while green up-quarks form Cooper pairs with red down-quarks, with the gap matrix A^. A priori, 
A^. As required by the overall antisymmetry of the spin-zero gap matrix P|, the right-hand side of Eq. Il2;-i|l is 
symmetric in color-flavor space. 

The regular quark self-energy was computed in Ref. 0|. To leading order, 

E*(A-).E-(AT.|l.„..l.(^) . (24) 

This corresponds to a wave function renormalization factor in the quark propagator. In the QCD ga p eq uation, it 
leads to subleading corrections which modify the prefactor of the color-superconducting gap parameter [2El2a|- Since 
we ultimately do not consider the QCD gap equation, but the one in a simpler point-Hke four-fermion coupHng model, 
we neglect the regular quark self-energy Ij24|l in the following. 

In order to proceed we compute the full inverse quark propagator Il21a|l with the gap matrix II23II , which is diagonal 
in the color-flavor basis lO , 

[G±]-^ =diag([G±:ri-AfGo^jAf, [G^tr - Af [G^^r ~ A^ cf^ A^ 

[G^y - Af Af , [Gr,r\ [G^tr) ^ (25) 

where 

[Gtir' ^rK^±^i{lo■ (26) 

At this stage, we have set the constant gluon background field A'^ — (color neutrality can be achieved by adjusting 
the color-chemical potential /ig), and we have also neglected the small up- and down-quark masses. The elements 
of the full inverse quark propagator Ij25|l have a simple physical interpretation. Consider, for instance, the red-up 
element [G^"]~^. The presence of the color-superconducting condensate Af (consisting of Cooper pairs of red up- 
and green down-quarks) modifies the propagation of red up-quarks, such that a red up-quark can be converted into a 
charge-conjugate green down-quark which continues to propagate and is then converted back into a red up-quark by 
the charge-conjugate condensate Aj~. 

With Eq. (jSHJ and setting = 0, Eq. Il21b|l reads in the color-flavor basis Q 

/ 

S±^^ 

„± _ 





V 

where 

„j.ud _ _ .± j-d „±(i« _ _ zpd . ± ±ti „±ud _ _ ^^u .± j-d „j.du _ _ ^d .± j-u , , 

^ rg — "-^0 r^l g ' ^ rg — ^0 r^2 g ' ^ gr ^ ^0 g ^2 "-^ r J ^ gr ~ ^0 g^l r ■ l-^oj 

Neglecting effects from the breaking of SU (3)c due to Cooper pair condensation (these effects are of sub-subleading 
order in the QCD gap equation ^SIJ), the gluon propagator can be taken to be diagonal in adjoint color, D'^'l{K — Q) = 



0\ 









0/ 



(23) 



o\ 






0/ 



(27) 
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5'^^D^^{K — Q). Inserting this in the gap equation H22a|l and performing the sum over adjoint colors we obtain the 
following four equations by identifying the non-trivial elements of the resulting color-flavor matrix: 



A+(if) 


Q 


1 —-i-ud / ^\ 






Q) , 


(29a) 




Q 








Q) , 


(29b) 




Q 








Q) , 


(29c) 


A+iK) 


^2+^(0)- 






Q) ■ 


(29d) 



Q 

In order to determine and A^, only two of these four equations are necessary. These two equations can be 
combined to 



4 , T 



A+(if)-f3A+(if) = -g^-Y^^^E+T,{Q)l''D,AK~Q), 



(30a) 
(30b) 



We now have to determine ^^"g and ■^^gr • The Dirac structure of the color-superconducting gap matrices Af^ is 
conveniently written in terms of energy-chirality projectors 



V!{k) = -(1 + C75) (1 + 6707 -k) , 



(31) 



where c = ± stands for the right/left-handed projection, and e = ± denotes the projection onto states of posi- 
tive/negative energy. With these projectors the gap matrices can be written as 



AtiK)=J2^ntiK)V!{k), A-iK)^J2^nriK)Vl!ik)^ n=l,2. 



(32) 



With the definition of the quasiparticle energy for positive and negative energy states (e = ±) 

elifi,(b)^^ik-efif + \<P\\ (33) 
we can write the components of the full quark propagator of the quark colors participating in Cooper pairing as 



where 



G+"(if) 


c,e 


G+t{K) 


-E 

c.e 


G+;{K) 


= E 

c.e 




-E 

c,e 



{ko 




(A,' 


hl)f 




^c(k) 






(fco 
















{ko 




(A,( 













■%{ko-5^i^f-[el (A,0i^)r 



[Go ;r'w , 

[GoT'iK) , 



5fii 

5^2 



2 2 
- 1 



Me , Ms 

M - ^ 



6 2v^ ' 



^ - § (M3 - Me) , 



2 

t4 - M^ 1. ^ . 

TT^ = 77 (M3 + Me) 



(34a) 
(34b) 
(34c) 
(34d) 

(35a) 

(35b) 
(35c) 
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Here we used Eq. ^ . Labelling the components of the full quark propagator with a color and flavor index is slightly 
misleading. For instance, red-up and green-down quasiparticles are both admixtures of red up and green down quarks. 
Both quasiparticles have the same Fermi surface fi, and only their dispersion relations ko — — (/i, (j)^^) ± Sfii differ 
by 2 

With Eqs. l|S3l and (EH and the relation G^A^G"^ = G^A^Gq, Eq. ^ becomes 



E 

c,e 

E 

c,e 

E 

c,e 

E 







Vc{±K) 


{ko 




[4 (A, 01^)]' 






^2t{±K) 


{ko 


T Sfi2)' - 






^|c(k)<? 


^2c(±i^) 


{ko 


± Sfi2)^ - 


[4 (A, 02c)]' 




^^c(k)<? 


^1?(±^) 


{ko 




[4 (A, 01^)]' 



(36a) 
(36b) 
(36c) 
(36d) 



Here we assumed that (j)j^'^{—K) = (j)^'^*{K). Taking the gluon interaction to be point-Hke, D^,y{K — Q) = — gf_i„/A?, 
A = const., the Dirac structure of the gap equations Ij3()ll can be projected out to yield gap equations for the gap 
functions 4>^'^^{K). It turns out that the gap equations for different energy, e = ±, and chirality, c = ±, projections 
decouple and have the same form. We therefore omit the indices e, c in the following. Due to our assumption of a 
point-Hke gluon interaction, the gap function is also independent of K. The gap equations l(3(Hl assume the simple 
form 



3 01 -I- 02 = 



3 A2 y 
3 A2 y 



EE 

Q e 

EE' 



{qo - Sfii 



[eq(M,0i 



1 2 ' 



Q e {qO+ 5^12) 

We now perform the Matsubara sums with the help of the relation 



T 



1 



The resulting gap equations are 

01+3 02 = 



1 



tanh 



5^1 



2T 



tanh 



- [eq(M,02)]' 
2T 



1 

2^ 



e{e~M) 



37r2 A2 ^ 



dq q^ 



02 = 



= -^|2E/ d'Z? 



eq(M,0i) 



eq(M,02) 



0(e^(M,0i)-|<5/xi|) 

^(6^(M,02)-|<5M2|) 



(37a) 
(37b) 

(38) 

(39a) 
(39b) 



Here we introduced a cutoff k to render the momentum integral finite. 

The equation of state for 2SC matter, i.e., the pressure as a function of temperature and chemical potential, results 
from the relation 



T 

P2SC = -j^ f 



(40) 



where F* is the value of the effective action JHJ at the stationary point determined by Eqs. ltTT|l . Since we do not 
consider the gluons as dynamical degrees of freedom, the first three terms in Eq. ijHJ can be omitted. The last two 
terms in Eq. ijHJ can be simplified with the help of the Dyson-Schwinger equation The final result reads 



P2SC 



1 T 

2 V 



Trln5"^ - - Tr (S 5) 



(41) 



where the propagator S obeys the Dyson-Schwinger equation Ijl5|l . Performing the trace over Nambu-Gor'kov space, 
the second term can be written as 



(42) 
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where the trace on the right-hand side runs over space-time, color, flavor, and Dirac indices. Since we neglected the 
regular self-energies in the derivation of the gap equation, to be consistent we also have to drop the first two terms 
in Eq. 1II2I1 . Fourier-transforming into momentum space and using the gap equations lf22ll , we obtain 



1 T 
4 V 



Tr (E S) 



^2 j,2 



^ Tr [rT,^-{K)^'^T^^+{Q)+^^TjE+{K)Yn^-{Q)\Df,{K^Q), (43) 

K,Q 



where the trace on the right-hand side runs only over color, flavor, and Dirac indices. We now insert the local, 
instantaneous gluon propagator D'^'l{K — Q) — — 5°"^ g^iu j Is? ^ and sum over a, h and ^, v. Then we perform the trace 
over color, flavor and Dirac space with the help of Eqs. lf27ll . Due to the point-like gluon interaction the sums 

over K und Q separate. These sums can be simplified with the help of the gap equations ll37ll . The final result is 



i I Tr (S 5) = ^ ^ -f 3 (0t02 -I- 010^)] 



The first term in Eq. l|lT|l is straightforwardly evaluated as 
1 T 



2 V 



■TrlnS'-i = 



T 
V 



I n=\a j=± 



y2 



In 



ckKN0)]2 



1 + exp - 



e^(Ai^0) + 2rin 
e^(M^,0) + 2Tln 



J^2 
T 

4(m^o) 



1 + exp 
1 exp 



T 
T 



In 



/r2 



J^2 



(44) 



(45) 



Subtracting the contribution from the vacuum and taking the limit T ^ leads to 



i-Trin5-i = i,y y 

2V 7r2 ^ ^ 



dk 



3 ek(M,0n) 



127r2 



(46) 



To obtain the pressure Ij41|l for color-superconducting quark matter with two flavors we have to add Eqs. Ij44|l and 

In Sec, mil we shall consider compact stellar objects which have to be neutral with respect to electric charge. In order 
to achieve this, we have to add the contribution of electrons to the pressure l|lTll of our two-flavor color superconductor. 



where we neglected the small electron mass. If the chemical potential for strange quarks, /if — fi + fie/i + fJ-3Tl^ + fJ'sT^i, 
exceeds the strange quark mass, nis, we also have to include strange quarks into our consideration. We assume them 
to be non-interacting, which leads to the following additional contribution to Eq. ()41ll . 

with the Fermi momentum kpt — (/if^ — m^y/^ and the energy — (fc^ + miy^^. Strange quarks also serve to 
neutralize the large positive electric charge of a system of up and down quarks. Consequently, we expect the electron 
density to be reduced once strange quarks are present in the system. The total pressure of our system is the sum of 
Eqs. gU, gZl, and M . 



P = P2SC + Pe+Ps - B 



(49) 
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Here, we have also subtracted the pressure of the perturbative vacuum in the form of the MIT bag constant B 
This will prove essential to obtain bound stars of finite radius. 

Compact stellar objects are not only neutral with respect to electric charge but also with respect to color charge. 
The neutrality conditions read 

ne = -^^0, (50a) 

713 = ^ - , (50b) 

ns = -^=0, (50c) 

where is the total electric charge density and 713, 713 are the color charge densities. It is straightforward to see that 
the solution of Eq. Ij50b|l is ^3 = 0. The reason is that in a two-flavor color superconductor, SU (3)c is broken to SU{2)c. 
One of the generators of this residual SU{2)c symmetry is T3 and, consequently, the associated color-chemical potential 
has to vanish (otherwise, SU{2)c would be broken explicitly). We therefore have /i3 = irrespective of whether we 
enforce color neutrality or not. >From Eqs. ll35|l we then conclude Sfj. = = —6^2 = — A*e/2. Inserting this result 
into the gap equations Ij^ijll we read off that the only possible solution is, 4> = ipi = —4>2- This greatly simplifles the 
gap equations; there is only a single gap equation for 

3^'7=E/^dqq2-^^0«(/2,</,)-|<SH) , (51) 

and the expression Il44|l becomes 

-i|Tr(E5) = -3^H^ (52) 
Finally, in order to solve the Tolman-Oppenheimer-Volkoff equation, we need the energy density, which at T = reads 

e = + /i3?l3 + ^8"8 + Me"e - P ■ (53) 

Here, n = dp/dfj, is the quark density. 

III. RESULTS 

In this section we numerically solve the gap equation Ij51|l and compute the equation of state II49|I . With Eqs. 
and ll53)l we then solve the TOV equation to obtain the mass-radius relation for quark stars. Unless mentioned 
otherwise, in order to compare our results to those of Ref. [2l|, we use the following values for the parameters of our 
model, 

= 0.1407 GeV , (54a) 
ffVA^ = 45.1467 GeV-2 , (54b) 
K = 0.6533 GeV , (54c) 
B^/" = 0.17 GcV . (54d) 

In Fig. [21 we show the total pressure lli9)l . normalized to its value for </) = = /ig = 0, as a function of the quark- 
chemical potential for various cases. The contribution of strange quarks is omitted, ps = 0, and the bag constant 
is set to zero, B = 0. One observes that the color-superconducting state with (p > has a larger pressure than the 
normal-conducting state with (/) = 0. Consequently, the color-superconducting state is energetically preferred. The 
difference in pressure between color-superconducting and normal-conducting states is proportional to the value of 
the gap (squared), cf. Fig. 31 The constraint of electric-charge neutrality reduces the pressure. (The effect on the 
pressure when imposing color-charge neutrality is negligibly small.) This constraint is necessary to obtain stable stars. 
Electrically charged stars would explode because of the repulsive Coulomb force. (One also has to impose color-charge 
neutrality because color-charged stars cannot exist due to confinement.) 

In Fig. the same cases are shown as in Fig. [2l now including the contribution of strange quarks. Again, the 
color-superconducting state is energetically preferred over the normal-conducting state. The difference to the previous 
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case is that, as soon as /i exceeds m^, strange quarks partially assume the role of electrons to ensure electric-charge 
neutrality of the system. Consequently, the electronic contribution to the pressure is reduced, cf. Fig. El and the 
reduction of the pressure when imposing electric-charge neutrality becomes smaller. For large values of the quark- 
chemical potential, ^ rris, the amount of electrons necessary to make the system electrically neutral becomes 
negligibly small. Therefore, the full and dashed lines, as well as the short-dashed and dotted lines in Fig. El approach 
each other. 

In Fig. 21 we show the gap as a function of the quark-chemical potential. For small values of the quark-chemical 
potential the gap vanishes. At larger values of /i, the gap increases, until the quark-chemical potential approaches 
the value of the cutoff k, where due to restricted phase space the gap starts to decrease. The conditions of electric 
and color neutrality Il5flajl and Il5flc|l tend to decrease the value of the gap. This decrease is especially pronounced 
in the case without strange quarks (dashed Hue), where the electro-chemical potential is large, see Fig. O Including 
strange quarks, the electro-chemical potential becomes smaller and, consequently, the gap becomes larger (dotted 
line) . Note the kink in the dotted line in Fig. 01 To the left of the kink, the system is in the so-called "gapless" color- 
superconducting phase discussed in detail in Ref. 0|, while to the right it is in the standard color-superconducting 
phase without gapless modes. For the dashed Hne, the system is always in the gapless color-superconducting phase. 

Figure shows the value of the electro-chemical potential /Ze as a function of /i as obtained from enforcing the 
neutrality conditions Ij50a|l and Ij50cll . For normal-conducting quark matter without strange quarks, electric-charge 
neutrality requires that the electron density increases proportional to the quark density. Consequently, fie is a lin- 
early rising function of fj, (full Hne). In the presence of strange quarks the amount of electrons required to achieve 
electric neutrality is smaller. Therefore, fie decreases as soon as /if > m^. For an electric- and color-charge neutral 
color superconductor either with or without strange quarks, the electro-chemical potential increases substantially as 
compared to a normal conductor. (Only with strange quarks and for very large values of fJ. close to k, the value of fie 
in the color-superconducting phase is smaller than in a normal conductor and may even become negative. We perceive 
the latter to be an artefact of approaching the limit of phase space.) 

In Fig. we show the value of the color-chemical potential fis as a function of fi for electric- and color-neutral 
color-superconducting quark matter with and without strange quarks. (Normal-conducting matter is automatically 
color neutral in the thermodynamic limit.) The values of fis necessary to make the system color neutral are much 
smaller than the values of fie required for electric neutrality. This can be understood in the weak-coupling limit where 
the gap is much smaller than the chemical potential, (j) ^ fi exp(— 1/q^) <^ fi. In this limit, the value of fis. required 
to achieve color neutrality is parametrically of order fi^ ^ 

With Eqs. and (IKHt we now solve the TOV equations 



The resulting mass-radius relations are shown in Fig. [3 We observe that the influence of color superconductivity on 
the mass-radius relation is at most on the order of a few percent, in agreement with the results of Ref. This was 
to be expected, since superconductivity is a Fermi-surface phenomenon, while the equation of state which determines 
the mass and radius is sensitive to the whole Fermi sea. To be more precise, the relative change in the pressure due to 
superconductivity is of the order (j)^ / fi^ . For (j) fi this is a tiny effect. Indeed, as can be seen from Fig.Q including 
strange quarks has a much larger effect on the masses and radii. The reason is that adding strange quarks to the 
system affects a relative change of the pressure by a term ~ [fig/ fi)'^ ~ C'(l). 

The observation that color superconductivity does not have an effect on the mass-radius relation of a quark star 
for 4> <^ fi immediately leads to the question how large the gap has to become in order to see an appreciable change 
in either the mass or the radius of the star. To our knowledge this question has so far not been addressed in the 
literature and is the main motivation for our current study. This question can be answered by artificially increasing 
the coupling constant g^/A^, such that the solution of the gap equation ll5T|l increases as well. Of course, such a 
modification will change the vacuum properties, such as the pion decay constant, within our NJL model. However, at 
nonzero baryon density, this model is in any case not very realistic as it neither describes saturation of nuclear matter 
in the ground state nor does it feature a hadronic phase where chiral symmetry is broken. Here, we consider the NJL 
model just as an effective model to describe quark-quark interactions at high baryon density. It is then certainly a 
permissible and interesting question to ask how the coupling strength affects the value of the color-superconducting 
gap parameter and quark star properties. 

In Fig. (HI we show the solution of the gap equation when decreasing the coupling constant by a factor of two (full 
line) as well as increasing it by a factor of 3/2 (short-dashed line) and a factor of two (dotted line), respectively. In 
the latter case we observe the interesting phenomenon that the gap is non- vanishing even when = 0. At first sight 



dp 
dr 



[p{r) + e(r)] [M (r) -I- 47rr3p(r)] 
r[r - 2M(r)] 



(55a) 




(55b) 
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this seems surprising but one can readily convince oneself that the gap equation Ij51|l indeed has a non-trivial solution 
for 11 = 0, if the coupling constant exceeds the critical value = Stt^/k^ ~ 69.37 GeV~^. The situation where the 

gap was non- vanishing in the vacuum was also analyzed in Ref. [sD]- 

In Fig. El we show the mass-radius relations calculated with the equations of state corresponding to the gaps shown 
in Fig. IHl As expected, a decrease of the coupling constant by a factor of two does not appreciably change the mass- 
radius relation, but multiplying the coupHng constant by a factor 3/2 already leads to an increase of maximum mass 
and radius by ^ 20%. For the coupHng constant which is a factor of two larger than the default value we find that 
the quark star also doubles in mass and radius. The reason is that the equation of state changes appreciably for gaps 
of order ~ 300 MeV, cf. Fig.|Hl since then 4> ^ ^. 



IV. CONCLUSIONS 



In this paper we have investigated color-superconducting quark matter in the so-called 2SC phase, i.e., a color 
superconductor consisting of massless up and down quarks, where quarks of red and green color form anti-blue Cooper 
pairs. Starting from the CJT formalism, we have derived the gap equation and the pressure. By adding electrons and/or 
strange quarks, we have also studied this system under the constraints of electric and color neutrality. For electric- 
and color-neutral systems, we have solved the TOV equation in order to determine the mass-radius relation for quark 
stars. We confirmed the result of Ref. |3| that color superconductivity does not substantially alter the mass and 
the radius of a quark star, if the coupling constant is chosen to reproduce vacuum properties such as the pion decay 
constant. The reason is that superconductivity has an effect of the equation of state which is proportional to 
which is small if <C ^. 

We then asked the question how large the color-superconducting gap parameter has to be in order to see an 
appreciable effect on the mass and the radius of a quark star. To this end, we artificially increased the coupling 
constant which has the effect of increasing the color-superconducting gap parameter. We found non-trivial solutions 
of the gap equation in the vacuum, i.e., for /i = 0, when the coupHng constant exceeds a critical value which is directly 
proportional to the cutoff parameter of our NJL-type model. For gaps of the order of 300 MeV the mass and radius 
of a color-superconducting quark star was found to be twice as large as for a normal-conducting quark star. While 
this is per se an interesting result, such quark stars are still of the same size and mass as ordinary neutron stars. It 
is thus impossible to decide whether a compact stellar object consists of normal-conducting or color-superconducting 
quark matter, or simply of ordinary neutron matter. 
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FIG. 2: The pressure as a function of the quark-chemical potential, normalized to its value for vanishing gap and fie = Us = 0. 
The contribution of strange quarks is omitted. The bag constant B is set to zero. Full and dashed lines are obtained for normal- 
conducting quark matter, = 0. Short-dashed and dotted lines are for color-superconducting quark matter, (p > 0. Full lines 
and short-dashed lines are computed without the constraints of electric- and color-charge neutrality, = Ms = 0. Dashed and 
dotted lines are computed for an electric- and color-charge neutral system, /le / 0, /is / 0. 
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FIG. 3: The same as in Fig. El but with the contribution of strange quarks. 
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FIG. 4: The solution <^ of the gap equation 15111 as a function of the quark-chemical potential jj.. The full line is for jj.e — ~ 0. 
For the dashed and dotted lines fie and fis are determined from the conditions Il5f)all and Il5flc|l of electric- and color-charge 
neutrality. The difference between dashed and dotted lines is that for the former the contribution of strange quarks is neglected, 
while for the latter it is included. 
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FIG. 5: The electro-chemical potential fie as a function of the quark-chemical potential jj. when enforcing the conditions of 
electric- and color-charge neutrality l|5f)a|l and l|50c|l . The full line and the dashed line are for normal-conducting quark matter, 
while the short-dashed and dotted line are for color-superconducting quark matter. The full line and the short-dashed line are 
obtained neglecting strange quarks, the dashed and dotted line are obtained including strange quarks. 
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FIG. 6: The color-chemical potential /is as a function of the quark-chemical potential /i, for electric- and color-neutral color- 
superconducting quark matter without strange quarks (full line) and with strange quarks (dashed line). 
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FIG. 7: The mass-radius relation for electric- and color-neutral quark stars. Full and dashed lines are for normal-conducting 
quark matter, short-dashed and dotted lines are for color-superconducting quark matter. Full and short-dashed lines are 
computed without, dashed and dotted lines are computed with strange quarks. 
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FIG. 8: The gap as a function of the quark-chemical potential /i. Full line: the coupling constant g^/A^ is reduced by a factor 



of two as compared to the default value Ij54hll . for which the corresponding gap is shown as the dashed line. Short-dashed line: 
the coupling constant is 3/2 times the default value; dotted line: the coupling constant is twice the default value. All gaps are 
for electric- and color-neutral matter including strange quarks. 
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FIG. 9: The mass-radius relation of electric- and color-neutral quark stars containing strange quarks. The full line is for normal- 
conducting quark stars, all other lines are for color-superconducting quark stars. Dashed line (almost indistinguishable from 
the full line): the coupling constant is 50% of its default value. Short-dashed line: the coupling constant assumes its default 
value. Dotted line: the coupling constant is multiplied by a factor 3/2. Dash-dotted line: the coupling constant is twice as large 
as the default value. 



